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SUMMARY 

In the walls of heat exchangers that are composed of noncircular 
passages, the temperature varies in the circumferential direction 
because of local variations of the heat-transfer coefficients. A pre- 
diction of the magnitude of this variation is. necessary in order to 
determine the region of highest temperature and in order to determine 
the admissible operating temperatures. 

A method for the determination of these temperature distributions 
and of the heat-transfer characteristics of a special type of heat 
exchanger is developed. The heat exchanger is composed of polygonal 
flow passages and the passage walls are heated uniformly by internal 
heat sources. The coolant flow within the passages is assumed to be 
turbulent. The circumferential variation of the local heat- transfer 
coefficients is estimated from flow measurements made by Nikuradse, 
postulating similarity between velocity and temperature fields. Cal- 
culations of temperature distributions based on these heat-transfer 
coefficients are carried out and results for heat exchangers with 
triangular and rectangular passages are presented. 


INTRODUCTION 

The conventional recuperative type of heat exchanger consists of 
passages for two liquids or gases separated by a heating surface. 

Heat from an outside source is carried with a fluid flowing through 
one of the passages and is transferred in the heat exchanger to a 
second fluid flowing through the other passage. Very often such a 
heat exchanger is composed of a large number of tubes, with the two 
liquids flowing inside and over the outside of the tubes, respectively 

The regenerative type of heat exchanger has passages for one 
fluid only. During the heating period, heat from an outside source is 
carried to the heat exchanger by a hot fluid and is stored within the 
solid walls of the passages. This heat is then given off to a cold 
fluid, which passes through the heat exchanger during the cooling 
period. 
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A heat exchanger is considered in this report that differs from the 
regenerative type only by the fact that the heat is generated by heat 
sources within the passage walls and is transferred to a coolant flowing 
continuously through the passages. The passage walls of this heat 
exchanger are assumed to be flat plates assembled to form a honeycomb. 
Thus the flow passages formed in this manner have a polygonal cross 
section. A cross-sectional view of a typical heat exchanger of this 
type is shown in figure 1. The exchanger is composed of a number of 
plates a, which form the coolant passages b. The flow of the coolant 
is normal to the cross section shown. High temperatures may be anti- 
cipated near the corners c of the passages, inasmuch as the rate of 
heat transfer there is expected to be poor. A theoretical investigation 
of the temperature distribution in such a heat exchanger was made at 
the NACA Lewis laboratory and is presented herein. 

The basis for this investigation is a knowledge of local heat- 
transfer coefficients in passages of noncircular cross section. Some 
information is available on the average heat- transfer coefficients in 
such tubes (references 1 and 2). The result of these investigations 
is essentially that the expressions derived for the heat-transfer 
coefficients in circular tubes apply for other cross sections as well, 
provided the diameter is replaced by the hydraulic diameter. (The 
hydraulic diameter is defined as four times the cross-sectional flow 
area divided by the circumference of the passage.) 

A knowledge of local heat-transfer coefficients in noncircular 
passages is important not only in the present problem, but also for 
several other engineering problems such as the determination of local 
temperatures in the walls of air-cooled turbine blades. No information 
on local heat- transfer coefficients was found, however, in the avail- 
able literature. These values are therefore estimated from flow 
measurements made by Nikuradse (reference 3) on the basis of the sim- 
ilarity between temperature and velocity fields. 


SYMBOLS 

The following symbols are used in this report: 

A cross-sectional flow area, (sq ft) 

A* ratio of wall area to flow area, (dimensionless) 

C internal circumference of passage, (ft) 

C* C/D, (dimensionless) 
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specific heat at constant pressure, Btu/(lb)(°F) 
hydraulic diameter, 4A/C, (ft) 
acceleration due to gravity, (ft)/(sec2) 
local heat- transfer coefficient, Btu/(sec)(sq ft)(°F) 

average heat-transfer 
Btu/(sec)(sq ft)(°F) 

h/ h, (dimensionless) 
thermal conductivity of wall material, Btu/(sec) (ft) (°F) 
thermal conductivity of coolant, Btu/(sec) (ft) (°F) 
k/kg, (dimensionless) 
residual value, (dimensionless) 

Nusselt number, .hD/kg, (dimensionless) 
coordinate normal to passage wall, (ft) 
n/D, (dimensionless) 

Prandtl number, v/a, (dimensionless) 
pressure, (lb)./(sqft) 

local rate of heat transfer, h0, Btu/(sec)(sq ft) 
h*0*, (dimensionless) 
radial coordinate, (ft) 

Reynolds number, • uD/v, (dimensionless) 

rate of internal heat generation, Btu/(sec)(cu ft) 

wall thickness, (ft) 

s/D, (dimensionless) 


coefficient, — h dx, 
C 1 ' 
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T local total temperature of coolant, °P 

Tg bulk total temperature of coolant, °F 

t local wall temperature, °F ' 

t* tk/rD^, (dimensionless) 

u,v velocity components in x- and y-direction, respectively. 

(ft/sec) 

x, y, z Cartesian coordinates 

x* x/d, (dimensionless) 

\ 

7* y/D, (dimensionless) 

z* z/D, (dimensionless) 

a thermal diffusivity, kg/pgCp, (sq ft)/(sec) 

3 angle subtended by two adjacent sides of polygonal passage, (deg) 

A increment of length, (ft) 

A* a/D, (dimensionless) 

6 , increment or difference 

turbulent diffusivity of heat, (sq ft)/(sec) 
turbulent diffusivity of momentum, (sq ft)/(sec) 

9 temperature difference, t-Tg, °F 

9* 9k/ rlfi , ( dimens ionle s s ) 

V kinematic viscosity, (sq ft)/(sec) 

p mass density, (lb) (sec^)/(ft)^ 

T w local wall shear stress, (lb)/(sq ft) 

T w average wall shear stress, ^ 
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Subscripts: 

c conditions for equivalent circular tube' 
m conditions at center of flow passage 

n conditions normal to passage wall 

s conditions at wall surface 

Superscript: 

* dimensionless quantity 


ASSUMPTIONS 

If this analysis were to be made without any simplifying assump- 
tions, the simultaneous solution of the equations of motion of the 
coolant and of the heat flow in the coolant and in the passage walls 
would be required. The following assumptions are made in order to make 
these equations amenable to solution without seriously curtailing the 
results of the analysis. 

(1) The Prandtl number of the coolant used in the heat exchanger 
is in the neighborhood of 1. This condition is well fulfilled by gases 
and by water above a temperature of 200°F, excluding the neighborhood 
of the critical point. 

(2) The passages of the heat exchanger are long enough so that in 

the cross section investigated the flow is fully developed, which, means 
that the velocity profile does not change its shape in the direction of 
the tube axis. , '• 

(3) The rate of heat generation in the walls of the heat exchanger 
is uniform. As a consequence of this condition, the temperature within 
the coolant and the walls increases linearly in a downstream direction 
provided the flow is thermally developed. For a fluid with a Prandtl 
number of 1, the points of thermal and velocity development in a tube 
practically coincide provided that the heating of the tube starts at 
the entrance section. 

(4) The temperature gradient along the tube axis is assumed small 
as compared with the gradients in any cross section of the passage. 
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(5) The thermal conductivity of the solid material is large as com- 
pared with the thermal conductivity of the coolant. This condition is 
always fulfilled for metal walls regardless of the type of coolant used 
provided that assumption (1) applies and for nonmetallic walls if the 
coolant is a gas. For the case of nonmetallic walls and liquid 
coolants, the applicability of the calculations presented in this 
report must be checked in each individual case. Furthermore, the tem- 
perature differences within the passage wall at any one cross section 
are postulated to be small as compared with the temperature differences 
between the wall and the core of the coolant. As a consequence of the 
assumption listed in this paragraph, the heat transport within the 
coolant, normal to the tube axis and parallel to the walls, is small 

as compared with the heat conduction within the walls and may therefore 
be neglected. The heat transport within the coolant, normal to the 
tube axis and normal to the walls, is of course taken into account. 

(6) The turbulent diffusivity of momentum € M and of heat Cg 
are equal. 


FLOW IN TUBES WITH NONCIRCULAR CROSS SECTION 

A thorough investigation of the flow through tubes with noncir- 
cular cross sections was made by Nikuradse (reference 3). This investi- . 
gation is used in the present report as the basis for estimating local 
heat-transfer coefficients. 

The flow of water through tubes of several different shapes, as 
shown in figure 2,' was investigated in reference 3. Three of the 
tubes had triangular cross sections; an equilateral triangle, an 
isosceles right triangle, and a right triangle with the sides enclosing 
the 90° angle, having a ratio of 1 to 2.32. One tube had a trapezoidal ' 
cross section and two tubes were circular with one and two grooves, 
respectively. In a second report Nikuradse investigated tubes with a 
rectangular cross section; however, only a summary of this report is 
available (reference 4). The hydraulic diameter of the passages varied 
from 0.3 to 0.6 inch and the length- to-hydraulic-diameter ratio varied 
from 100 to 200. The Reynolds number based on the hydraulic diameter 
and the mean velocity ranged from 77,000 to 120,000. Because it is 
known that the shape of a turbulent velocity profile changes only 
slightly with Reynolds number, the results of the calculations should 
be applicable for a fairly large range of Reynolds numbers in the 
turbulent region. The velocity profiles were measured with a small 
total-head tube in a cross section near the downstream end of the tube. 
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Results of reference 3 for the isosceles right triangle are shown 
in figures 3 to 5. Similar results for all -other cross sections can 
be found in the original report. In figure 3 the lines of constant 
velocity are presented as contour lines. These contours are indented 
near the center of each side of the passage. The conclusion of refer- 
ence 3 is that these indentations indicate a secondary flow normal to 
the main flow direction. A qualitative sketch of this secondary flow 
is presented in figure 4. Prandtl (reference 5) has ascribed the sec- 
ondary flow to a turbulent mixing motion within the fluid, which is 
more intense in the direction parallel to the wall than normal to it. 
This secondary flow tends to equalize the velocities and temperatures 
within any cross section of a noncircular tube and is therefore favor- . 
able for the similarity consideration that will be used to deduce the 
temperature field from the measured velocity field. On the other 
hand, because no quantitative knowledge of the secondary flow exists, 
an exact theoretical calculation of the temperature field in noncir- 
cular passages is impossible. The local wall shear stresses were 
computed by Nikuradse, using the measured velocity profiles together 
with the assumption that the Blasius pipe resistance law of the tur- 
bulent velocity profile for circular tubes applies also for noncir- 
cular passages on normals to the walls (fig. 5). A check by Nikuradse 
of the calculated average shear stresses against the measured pres- 
sure drop showed good agreement. 


SIMILARITY BETWEEN VELOCITY AND TEMPERATURE 
PROFILES FOR Pr = 1 


It is well known that the velocity and temperature profiles in 
the boundary layer of a fluid with a Prandtl number of 1 are similar 
in shape in the absence of a pressure gradient. This similarity is 
immediately apparent from a comparison of the momentum and energy 
equations of the boundary layer for two-dimensional steady flow along 
a flat plate (reference 6): 


p ( u s + T If) = ^ [ pfc+e M ) ijl] 

p ( u m - t n) • 4 [ p(,,+£b) n] 


a) 

(a) 


where T represents the total temperature, the Prandtl number is equal 
to 1, and the specific heat is constant. The effect of internal fric- 
tion on the temperature profile is taken into account by basing equa- 
tion (2) on the total temperature. 
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Because it is assumed that the turbulent diffusivity of momentum 
€ m and of heat are equal', equations (l) and (2) are similar. 

The solutions of the equations, namely the velocity and temperature 
profiles, are therefore also similar, provided the boundary conditions 
are similar. The boundary conditions for the velocity field on a flat 
plate are: the velocity is 0 along the wall and has a constant value 

outside the boundary layer. Similar boundary conditions for the temper- 
ature field are: The temperature is constant both along the wall and 
in the main stream. 

It is therefore evident that similarity between velocity and tem- 
perature profiles in the boundary -layer over a flat plate exists in an 
exact mathematical sense. The same is not true, however, for fully 
developed pipe flow. The balance of the forces and of the heat energy 
on a stationary annular volume element with a radius E, thickness dR, 
and length (in the direction of the tube axis) dz, leads to the fol- 
lowing equations for fully developed flow in a circular tube: 

' 5 m [rp(v+€ m ) || = H (3) 

1 J? Jrp(cw-%) H] - P u § (4) 

where the heat-conduction term in the z-direction is neglected in equa- 
tion (4). It is apparent that equations (3) and (4) are not similar, 
even if the viscosity and diffusivity terms in the two equations are 
equal. The following type of analysis can be made, however: 

The turbulent diffusivity of momentum can be calculated from 

equation (3), provided that the velocity profile u = f(E) and the 
pressure drop dp/dz are experimentally known. If it is again assumed 
that €j| and €g are equal, the temperature profile can then be cal- 
culated from equation (4). Such an analysis, based on constant property 
, values and with internal friction neglected was carried out by Latzko 
(reference 7). The results of this analysis- were used to calculate the 
relation between the temperatures and the velocities for hydrodynamically 
and thermally developed flow of a fluid with Pr = 1. Figure 6 rep- 
resents a plot of -the local temperature -difference ratio 
(T-t s )/(T m -t g ) against the local velocity ratio u/u^j. It can be 

seen that this relation is very nearly linear. This linearity means 
that the velocity and temperature profiles are similar as far as 
practicality is concerned, even though they are not similar in a 
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strict mathematical' sense. Most of the recent theoretical investiga- 
tions of turbulent heat transfer are therefore based on the assumption 
of similarity between these profiles for a fluid with a Prandtl number 
of 1. 


It is known that the fully developed turbulent velocity profile 
in a circular tube is such that the essential part of the velocity 
change from 0 at the wall to a maximum value at the center occurs in 
a narrow strip around the periphery of the tube. The investigations 
of references 3 and 4 show that this condition is also true for tubes 
of polygonal cross section, as long as the included angles between two 
adjacent passage walls are not very small. Nikuradse further shows that 
the law for the velocity variation on a normal to the wall established 
for turbulent flow in circular tubes holds also for the noncircular 
passages in this region of essential velocity variation. It can there- 
fore be concluded that the velocity and temperature fields for fluids 
with a Prandtl number of 1 flowing through tubes of polygonal cross 
section are practically similar. 

An immediate result of the similarity between velocity and tem- 
perature profiles is the fact that the local wall shear stress deter- 
mined by the velocity gradient at the wall is proportional to the local 
heat-transfer coefficient determined by the temperature gradient at the 
wall. (The same result can be obtained from Reynolds’ analogy). The 
proportionality of these values, together with the knowledge that the 
average heat-transfer coefficients for circular and noncircular tubes 
with the same hydraulic diameter are equai, can be used to obtain local 
heat-transfer coefficients in noncircular tubes from the wall-shear- 
stress data of references 3 and 4. A generalization of this result 
is possible because, for a given Reynolds number, the shear stress, 
distribution in the vicinity of a corner of the passage depends pri- 
o marily on the conditions near that corner and not on the shape of 
the passage. 

The local heat-transfer coefficients, as presented in figure 7, 
were obtained by correlating the results of references 3 and 4. In 
this figure the ratio h* of the local heat-transfer coefficient to 
the mean value is plotted against the dimensionless distance from the 
corner x*. The included angle of two adjacent passage walls is the 
parameter for the curves. Some of these comer angles were represented 
on more than one of the cross sections investigated in reference 3, and 
for these angles the curves of h* against x* agree reasonably well. 

The application of figure 7 can best be explained with the aid of 
an example. Suppose it is desired to find the distribution of h* 
over the short leg of an isosceles right triangle, as shown in figure 8. 
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The two ends of the curve can immediately he transposed from figure 7 . 
The remainder of the curve is then extrapolated so that the area under 
the curve divided hy length AB is equal to unity. There is little 
freedom in the extrapolation of the curves, and the error produced by 
this extrapolation is of the order of 5 percent. 


DIMENSIONLESS VARIABLES 


With the distribution of the local heat-transfer coefficient around 
the circumference of the passages established, the only problem remain- 
ing to be solved is the calculation of the heat- conduction process 
within the passage walls. Before this problem is taken up, an investi- 
gation to determine on what dimensionless moduli the temperature dis- 
tribution in a heat exchanger of the type under consideration depends 
is made. This investigation may be useful as a basis for experimental 
investigations. None of the simplifying assumptions, as summarized 
previously, is necessary for the development in this section of the 
report . 


The heat conduction within the solid walls of a heat exchanger 
with internal heat generation is described by Poisson's equation 




b 2 t d 2 t d 2 t\ 

,3x 2 + dy 2 + Sz 2 / + r = 


(5) 


The number of parameters for this equation can be reduced by the use of 
the following dimensionless values: 





( 6 ) 


With these values equation (5) can be rewritten 


d?t* + ^E + ^ + l = 0 

cbc * 2 Sy * 2 dz * 2 


(?) 


The temperature field that results under the influence of the internal 
heat generation in the solid walls depends on how the heat is transferred 
from the wall surfaces to the coolant. The heat must be conducted within 
the solid material to the surface and from there into the coolant? Along 
the surface, therefore, the following boundary condition exists: 



( 8 ) 
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where n represents the direction normal to the surface. In terms of 
dimensionless variables this equation "becomes: 

*•©. ■ (s). . 

The heat transfer within' the. coolant "by conduction and convection 
is governed by an energy equation similar to equation (4) and the flow 
of the coolant is determined by a momentum equation (similar to equa- 
tion (3)) and the corresponding continuity equation. Textbooks on 
heat transfer (for example, reference 8) show that the dimensionless 
temperature described by the energy equation of the coolant depends, 
for low-velocity flow and constant property values, on the Reynolds 
number Re and the Prandtl number Pr. At high velocities and vari- 
able property values, there is an additional influence of the Mach num- 
ber and other dimensionless expressions characterizing the temperature 
dependency on the property values. By neglecting the last-mentioned 
influences and summarizing all the factors that influence the heat 
flow in the solid material and in the coolant, a functional relation 
of the following type can be deduced: , 

6 * = f (Re,Pr,k*,x*,y*,z*) . (10) 

where the temperature difference 9 * is introduced because only tem- 
perature differentials appear in the equations. The function f as 
expressed by this equation depends on the geometric configuration of 
the heat- exchanger passages. 

When it can be assumed that the heat-transfer process from the 
walls to the coolant is not influenced by the temperature distribution 
within the wall, then the number of factors influencing the problem 
can be considerably reduced. TJiis assumption is made throughout the 
calculations in this report and seems reasonable as long as previously 
mentioned assumption (5) holds. In this case equation (8) can be 
replaced by ( 

k (!s) “ - h (V T B> (11) 

' 'S 

where h is a known function of the flow parameters. By introducing 
the ratio h* of the local heat-transfer coefficient to the value 
averaged over the circumference of the passage, the Nusselt number 
Nu = hD/kg based on this average value, and by changing to the tem- 
perature difference 6, equation (11) can be transformed to 
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The dimensionless temperature vithin the solid walls of the heat 
exchanger is determined by equations (7) and (12) and' can be presented 
as a function of the following kind: 

6* = f (Nu g,x*,y*,z*) (13) 

This dimensionless temperature difference depends only on the dimen- 

sionless local coordinates and on the parameter Nu — • Temperature 

distributions and heat-transfer characteristics for any given geometric 
configuration are therefore presented as a one-parameter family of. 
curves. 


CALCULATION OF TEMPERATURE DISTRIBUTIONS WITHIN PASSAGE WALLS 

The problem of calculating temperature distributions within the 
passage walls can be classified in two general categories, depending 
on the relative thickness of these walls. If the dimensionless wall . 
thickness s* (see fig. 8) is small as compared with the length of a 
wall AB, then the temperature differences in the direction normal to 
the Wall surface are small as compared with the temperature difference 
in the direction x* parallel to the surface. Only the temperature 
differences in this parallel direction need therefore be considered. 
Hereinafter this special case of the problem is referred to as the 
"one -dimensional problem." On the other hand, if s* becomes large, 
then the heat flow in both the x* and s* directions must be con- 
sidered. This second and more general problem is referred to as the 
"two-dimensional problem." 

One-dimensional problem. - Consider a thin plate of thickness s 
that is a section of the heat- exchanger passage walls. Heat is gen- 
erated uniformly throughout the plate at a rate r. The local plate 
(or wall) temperature t is assumed to be varying in the x-direction 
only. The bulk temperature of the coolant is Tg and the thermal 
conductivities of the wall and coolant are k and kg, respectively. 

The local coefficients of heat transfer to the coolant above and below 
the plate are denoted by h and h', as shown in the following sketch: 
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The heat balance for an element of volume vith the dimensions dx 
and s and of unit depth is 


where 


Q+Qi = Q'+Qa+V 


^ dt 

Q = -kS ^ 


Q* = -ks 


d^t 
+ dx 

dx 2 




Qi = rs dx 

Q a + Qa’ = (h+h’)(t-T B ) dx 

With the preceding values of Q, equation (14) becomes 

- ^|1 (t-T*) + ? = 0 


dx 


B' 


(14) 


(15) 


Equation (15) can be expressed in terms of dimensionless variables, 
d 2 0* h*+h** 


dx*2 


s*k* 


Nu 0* + 1 = 0 


(16) 


This equation expresses the dimensionless wall temperature 9 * as 
a function of the dimensionless coordinate x* and a single parameter 
(h*+h'*)Wu/s*k* . For any given geometric configuration (which deter- 
mines h* and h'* as functions of x*), the dimensionless wall tem- 
perature is a function only of the average Nusselt number Nu, the 
dimensionless wall thickness s*, and the ratio of the thermal con- 
ductivities k* . 

In general, equation (16) must be solved by numerical means inas- 
much as h* and h'* are experimentally determined functions of x* . 
Several numerical methods of solution can be applied and two of these 
are considered here. The relaxation method (see, for example, refer- 
ence 2, pp. 365-379, or reference 9) has the advantage that it is easy 
to apply and that computational errors are immediately apparent. A 
second method of solution, based on the Runge-Kutta method (refer- 
ence 10) is presented in the appendix. ‘This method is preferable when 
results of high accuracy are desired. It is less convenient, however, 
than the relaxation method because computational errors are much more 
difficult to detect. 
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For the purpose of solving equation (16) by the relaxation method, 
the equation is first expressed in finite-difference form. Consider 
a grid, or net of points, placed into the wall, any two adjacent 
points being separated by a small but finite distance A*, as indi- 
cated in the following sketch: 



I— A* -j- A* 

The second derivative of 9* at an arbitrary point 0 can be 
expressed in terms of the temperatures 9* at this point and the two 
adjacent points as follows: 



’ ^2 ( e l* +0 2* - 2S 0*) (17) 

With this value for the second derivative, equation (16) becomes 

®1* + d 2* - e 0 * [ 2 + Jgr (fc*+fc’*)oA* 2 ] + A* 2 = 0 (18) 

The relaxation method, used for solving this equation, is subsequently 
discussed. 

Two-dimensional problem. - Consider a slab of homogeneous solid 
material with unit thickness, . as shown in the following diagram: 

Heat is again generated within the 
material at a uniform rate r. The 
heat balance for an element of vol- 
ume with the dimensions dx and 
dy and of unit depth is, 

V v Qi = v+v (19) 


1- dx-l 



x 
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A constant flow of heat through the slab normal to the plane of 
the preceding sketch does not influence this heat balance. Local devi- 
ations from this constant flow of heat are neglected. 

The individual terms in equation (19) are 


- - k s 


Q ' = -k dy 


(|S + aft ax) 


V 


Qy = -k to || 

. . /dt d 2 t . \ 

- kdI W + 55 a! 7 


'Qj. = r dx dy 

Equation (19) can therefore be written as follows: 

d 2 t + d 2 t _ r 
dx 2 dy 2 k 

Or, in terms of the dimensionless variables, 

ci 2 e* a 2 e» _ 

dx*** dy*^ 


( 20 ) 


( 21 ) 


If the slab is bounded by a coolant whose bulk temperature is Tg and 
if the local surface heat- transfer coefficient at a point s on the 
surface is h, then the boundary condition at that point is 



( 22 ) 


where n -is. the direction normal to the surface. Equation (22) can 
again be written in terms of the dimensionless variables 




Nu ha- 
le* 



( 23 ) 
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Equation (21) together with boundary condition (23) fully describes the 
two-dimensional problem. All physical variables are again grouped into 
a single parameter h* Nu/k*, which appears in the boundary condition. 
In order to apply the relaxation method of solution to these expres- 
sions, they are first converted to difference equations. 


For this purpose a rectangular grid is placed into the slab, as 
indicated in the following diagram: 

Adjacent net points are sep- 
arated by a distance A*. 

The derivatives of 9* at an 
arbitrary point 0 can be 

A* | I expressed in terms of the tem- 

perature function at surround- 
ing points. Thus the first 
derivatives are 


y* 


o 4 

—1 

1 

1 

M 

o 

o 

o 

o 

CM 

1 

o 2 

1 

_J 


I- A* 4 

(be*\ 

M 

(be*\ 

2 W */0 


X* 



(be*) . 

V-*: 

A* 

o\^/3 ' 

A* 

ft * p * 

(be*) . 

0 O *- 0 - 

A* 

oW 

A* 


and the second derivatives become 


d 2 ft* _ e l* + 0 5*~ 2e Q* 
dx* 2 A* 2 

d 2 0* _ 0 2* +0 4*~ 20 O* 

by* 2 . A* 2 • 

With these values, equation (21) becomes 

0 l* +0 2* +0 3* +0 4*- 40 O* +A * 2 = 0 - (24) 

The boundary condition can be evaluated by assuming that points 1, 
0, and 3 lie on the surface of 'the wall. Point 4 then lies in the 
stream and its temperature must be expressed by the normal derivative 
at point 0. 
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/<#*\ s 4 *-e 0 * 

4 W 0 A * 

or^ with the use of* the value of S0*/dy* as given by equation (23), 

V-V^-TT") (25) 

Along the boundary, therefore, the following equation applies: 

0 l* +0 2* +0 3*" 0 O* ( 3 + A *) + < A *) 2 = 0 (26) 

Equation (24) together with boundary condition (26) can again be 
solved by means of the relaxation method. 

Solution by relaxation method. - A heat exchanger composed of a 
large number of rectangular passages will be discussed. The walls of 
the heat exchanger are made of a homogeneous material in which heat is 
generated at a uniform rate r. Figure 9 represents a cut through the 
heat exchanger so that the flow of coolant is in a direction normal to 
the cut. In the discussion that follows, it is assumed that the geo- 
metry of the configuration is given in the dimensionless system of 
coordinates . 

A complete discussion of the relaxation method is not presented 
herein, inasmuch as it is generally available elsewhere. (See, for 
example, reference 9.) The essential features of the method can be 
outlined as follows: 

Suppose it is desired to solve a given finite-difference equation 
over a certain area of integration. The equation is of the following 
type: 

f(0*) =0 (29) 

The solution of equation (29) must also satisfy prescribed conditions 
at the boundary of the area of integration. 

First, it is necessary to select a number of net. points covering 
the entire area of integration. The distance A* between net points 
is arbitrary, with the accuracy of the final solution increasing as the 
distance between points is decreased. Next, values of 9 * are ass ume d 
at each net point. If by chance these assumed values of 9 * are the 
correct values, then they satisfy the appropriate finite-difference 
equation at all net points. 
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In general, however, the assumed values of the function do not 
satisfy the difference equation and the left-hand side of equation (29) 
is equal to some residual value N' instead of zero. At any given net 
point it is then nece.ssary to adjust 9 * in order to make, N vanish 
at that point. This adjustment of 9 * then also changes the residuals 
at adjacent net points. However, if this process of adjustment is 
started at the point at which the absolute value of N is greatest, 
and is then repeated for points at which the value of the residual is 
successively less, the correct values of 9* for the entire net 
eventually are obtained. 


In applying this method to the heat exchanger under consideration, 
it is first assumed that the wall thickness s* is small, so that the 
one -dimensional solution applies. 


The actual configuration to be discussed is shown in figure 10. 
Because points A, B, and C are points of symmetry, only the region 
bounded by these points need be considered. Furthermore, the distri- 
bution of temperature is also symmetric about point 0 with respect to 
AO and 0C. An enlarged view of the section under consideration is 
represented in figure 11. 

The section has been subdivided into a number of net points a, 
b, c, . . . h, each separated by a distance A* from adjacent points. 
The first step in the analysis is the assumption of temperatures at 
all net points, which can be done by setting up a balance between the 
total heat produced within the walls and the total outflow of heat 
from the walls. t In terms of the dimensionless variables, this. heat 
balance becomes 

^C* 

e*(h*+h’*) dx* = 1 (30) 

0 


JL_ Nu. 
C* k*s* 


If, for the initial assumption, h* = h** = l and 9 * is constant, the 
following expression is obtained: 


0 * = 


s*k* 
2 Nu 


(31) 

\ 


The value of 9 * given by equation (31) is assumed to exist at all 
net points. The appropriate difference equation at points b, c, d, f, 
and g is (from equation (18)) 

0 l* +0 2* _0 O* [ 2 + (h*+h’*) 0 A* 2 ] + A* 2 = N 0 (32) 
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where the subscript 0 refers to the point at which the equation is to 
be applied and subscripts 1 and 2 refer to adjacent points. At 
point a, which is influenced by three points, the following equation is 
required: 


20b*+ 8 f*- e a* Q + i*Ef (h»+f») a 4.2] + 4.2 = H a (33) 

At point e, which is a point of symmetry, the following expression 
applies: 

2e d*- e e* [ 2 + ^ (*r*+h'*) e A*2] + A*2 = N e (34) 

A similar expression applies at point h. With these equations, the 
residuals N are calculated at each net point. The point at which N 
has the largest absolute value is then selected and 0* at that point 
is adjusted so that the residual vanishes. The effect of this adjust- 
ment on the residuals at adjacent points is calculated with the aid of 
the appropriate finite-difference equation. 

The process is repeated for the point at which the next largest 
value of the residual appears. Eventually, if the process is repeated 
often enough, the residuals at all net points approach zero. The final 
adjusted values of 0* then satisfy the appropriate equations at all 
net points . 

These values of 0* should also satisfy the heat balance as given 
by equation (30). This equation can therefore be used to check the 
validity of the final temperature distribution. If the check is unsat- 
isfactory, it is necessary to continue the solution of the finite- 
difference equations by using a smaller net spacing. 

The same heat -exchanger configuration is now investigated without 
making the assumption that the wall thickness is small. The two- 
dimensional equations are therefore applied. An enlarged view of the 
shaded portion of figure 9 is shown in figure 12. The lines CD, FG, 
and GA are lines of symmetry, and there is no flow of hea t across 
these lines. The temperature distribution along line DEF will be 
symmetric about point E, but there will be a flow of heat across this 
line. A large number of net points is selected to cover the entire 
section of the configuration shown in figure 12. (Only a few of these 
points are presented in the fig.) An initial value for the temperature 
along lines AB and BC is obtained by setting up a balance between the 
total heat produced within the section and the flow of heat to the 
gas. In terms of the dimensionless variables this heat balance becomes 
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_ 1 _ 

C* 


'■'C* 

h*e* ax* 

J 0 


A* k* 
4 Nu 


(35) 


If, for the initial assumption, h* = 1 and 0* 
then the following expression for 6 * is obtained: 


is constant. 


0 * = 


& & 
4 Nu 


(36) 


It can now be assumed that the temperature as given by equa- 
tion (36) exists along the surface of the wall and that a somewhat 
higher temperature exists at internal points. The ass ume d temperatures 
are then adjusted by means of the relaxation method. 


The appropriate finite-difference equation for internal points 
such as point e, is written as follows: 

0 b* +e f* +e h* +0 d*- 40 e +A * 2 = N e (37) 

For points along the boundary, such as point b, the equation is 

e a *+ 0 e* +0 c*- 0 b* ( 3 + P h b* A *) + A* 2 « N b (38) 

and for points, on a line of symmetry, such as point f , 

0 c *+20 e *+e i *-4e f *+A* 2 = N f (39) 

The appropriate equation for points along KEF can be determined by 
taking advantage of the antisymmetry about this line. At point k for 
instance, the equation becomes 

e h*+ 0 1 *+ 0 n # + 0 j*-40k* + ^ 2 = N k (40) 

Expressions of the type (37) to (40) apply at all net points. The 
method of determining the actual temperature distribution is the same 
as the method outlined in the previous section. 

The final' temperatures along the surface of the wall can be 
checked with the aid of equation (35). If the check is unsatisfactory, 
a finer net spacing is required. 
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RESULTS 

One-dimensional solution. - The one -dimensional solution was . 
applied to heat exchangers composed of rectangular and triangular pas- 
sages, respectively. In each case the passages were staggered in order to 
minimize the expected hot spots. 

The rectangular configuration is represented by figure 10. The 
height-to-width ratio of each passage in this configuration is 1 to 5. 
Figure 1 represents the triangular configuration. Each passage in 
this configuration is an isosceles right triangle. 

It has previously been stated that, for each configuration, the 
temperature distribution in the passage walls is a function only of 
the parameter Nu/ s*k* . Accordingly temperature distributions were 
calculated for several values of this parameter, ranging from 0.1 to 
100. It is believed that this range of values includes all values 
actually encountered. 

Temperature distributions were calculated according to the method 
that was outlined in the previous section of this report and were then 
checked by equation (30). For this purpose local rates of heat transfer - 
q* = h * 6 * were calculated. According to equation (30), the mean value 
of q* when multiplied by 2Nu/k*s* should equal 1. Because the 
finite-difference method of solution is essentially an approximate 
method, the 'results were not expected to be exact. All results pre- 
sented in this section of the report, however, were held to an error 
of less than 6 percent. The final temperatures were multiplied by a 
constant scale factor in order to satisfy equation (30) exactly. 

Temperature distributions and local rates of heat transfer are 
presented for the triangular heat exchanger in figures 13 and 14 and 
for the rectangular heat exchanger in figures 15 and 16. In the tem- 
perature plots (figs. 13 and 15), the temperature difference ( t-Tg) 
is referred to the temperature difference (t-Tg) c for a thin- walled 

circular tube with the same hydraulic diameter and an internally 
heated wall with the same physical properties. The temperature dif- 
ference for the circular tube is given in the nondimens ional form by 
equation (31). Similarly, the ordinate in figures 14 and 16 refers 
the rate of heat transfer q to the rate of heat transfer from the 
wall of a thin circular tube q c . It is believed that the interpreta- 
tion of the temperature and heat-transfer curves is simplified by this 
final change of ordinates. Temperature distributions for the trian- 
gular configuration are shown in figure 13. These curves are presented 
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for values of Nu/s*k* of 0.1, 1.0, 10, and 100. As vas expected, the 
temperatures peak near the corners of the passage. The largest tem- 
perature differences are encountered for large values of the parameter 
Nu/s*k*. A large value of this parameter indicates a low conductivity 
of the wall material or a small wall thickness. Both of these condi- 
tions are conducive to low rates of heat exchange, within the wall. 

Local rates of heat transfer corresponding to the temperature distri- 
butions just discussed are presented in figure 14. 

Temperature and heat-transfer curves for the rectangular config- 
uration are shown in figures 15 and 16. The results for this config- ' 
uration are similar to those for the triangular configuration. 

In order to determine when the one -dimensional solution can be 
used, it is necessary to obtain some information on the temperature 
difference that exists across the passage walls on normals to the 
surfaces and to compare this temperature difference with the temper- 
ature differences along the surfaces. An estimate of the tempera- 
ture difference . across the wall may be obtained by calculating 

this value for a flat plate. The result of this calculation is 


5t 


n 



(41) 


The equation can be changed to the dimensionless values and yields 




(42) 


This temperature difference can again be referred to the temperature 
difference (t-Tg) c 


_ s*Nu 
(t-T B )c 4k* 


(43) 


The one-dimensional solution applies as long as this value is small as 
compared with the temperature differences presented in figures 13 and 15 . • 

Two-dimensional solution. - The two-dimensional solution for the 
temperature field within the heat exchanger depends on two parameters; 
namely, the dimensionless wall thickness s* and the value Nu h*/k*. 

In addition, the time required to obtain the solution for a special 
case is much longer than for the one-dimensional solution. Only one 
example was therefore calculated; namely, the temperature distribution 
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within the walls of a heat exchanger composed of rectangular passages, 
as shown in figure 9. 'The ratio of the two side lengths of the rec- 
tangle is 1 to 5. The ratio of the hydraulic diameter to the short side 
of this passage is 1.67 and the dimensionless wall thickness is 0.6. 
Figure 17 presents the results of the calculation using the relaxation 
method with a network of 88 points.- Lines of constant temperature 
(isotherms) are shown within the portion of the heat- exchanger walls 
that is shaded in figure 9. It may be seen that the heat flow within 
this wall is mainly in the direction normal to the wall surfaces. The 
maximum temperature differences on any normal to the surface are not 
very different from the value in a flat plate, as calculated in equa- 
tion (42). The temperature differences along the surface of the wall 
are appreciably smaller -than the temperature differences across the 
wall. It cannot therefore be expected that the values of the one- 
dimensional solution apply to this case. Actually, the temperature 
difference found in figure 17 along the surface of the wall" is con- 
siderably greater than for the corresponding one-dimensional case. In 
addition,, heat is generated within the corner area of the rectangular- - 
wall configuration and has to be conducted away along comparatively 
long paths. An additional temperature increase can therefore be found 
within this corner area. 


CONCLUSIONS 

A method to calculate the temperature distribution in a heat 
exchanger composed of noncircular flow passages with internally heated 
walls has been presented. 

Local heat-transfer coefficients along the circumference of the 
heat-exchanger passages were obtained from flow measurements made by 
Nikuradse, assuming similarity between the velocity and temperature 
fields. The heat-transfer coefficients, as determined in this manner, 
decrease sharply near the corners of the passages and vanish at the 
corners. This- decrease becomes more pronounced as the corner angle is 
increased. Near the corners these coefficients are essentially influ- 
enced only by the magnitude of the included angle of the corner . 

It was shown that the dimensionless temperature distribution 
within the passage walls depends on a single parameter, provided the 
dimensionless wall thickness is small. Numerical evaluations for 
triangular- and rectangular-passage configurations for a wide range 
of the aforementioned parameter were carried out. The results of 
these evaluations are presented so that temperature differences 
arising in the walls of heat exchangers with the investigated passage 
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shapes for any condition within the range of practical interest can he 
read off one of several curves . These curves show a temperature 
increase near the corners of the passages. This increase becomes more 
pronounced for high Nusselt numbers based on the average heat-transfer 
coefficient, for low wall- thickness ratios, and for low ratios of the 
conductivity of the wall to the conductivity of the coolant. 

The dimensionless temperature distribution for thick-walled heat 
exchangers depends on the wall-thickness ratio in addition to the 
property parameter. Because the time required to calculate temperature 
distributions in the whole field of interest determined by the two 
parameters is prohibitive, only a specific example was evaluated 
numerically. The method of calculation is presented in great detail, 
however, so that evaluations for other interesting cases can be 
carried out. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, October 4, 1950. 
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APPENDIX - SOLUTION BY EUNGE-KUTTA METHOD 
Equation (16) is of the following form: 

-f(x*)6* + 1 = 0 

dx*^ 

The boundary conditions are 



where a and b are points of symmetry. 


(Al) 


(A2) 


In general, a solution of equation (Al) could be obtained by 
assuming an initial value of 0* at point a in addition to the first 
of the boundary conditions (A2) and by working the solution towards 
point b using a numerical method of integration. If the boundary con- 
dition at point b is not satisfied, the solution has to be repeated 
with a new initial value. This trial-and- error process is tedious and 
can be avoided by splitting equation (Al) in such a manner as to keep 
an undetermined constant c in the solution. This constant is 
finally determined by satisfying the second boundary condition. 

Let 

e* « Sj* + c e 2 * (A3) 

d 2 e * 

-f(x*) ©i* + 1 = 0 

dx* 

-f(x*) e 2 * = 0 
dx* 



where 
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These expressions still satisfy equation (Al) and the first of the 
boundary conditions (A2). After, the solutions for 0-j* and 0g* have 

been numerically obtained, the constant c is determined so that the 
boundary condition at point b is fulfilled 



The final temperature distribution is given by equation (A3). 

The temperature distribution in a thin-walled heat exchanger with 
triangular passages was calculated in this manner (for Nu/s*k* = 10) 
and was compared with the corresponding curve in figure 13. The Runge- 
Kutta method (reference 10) was used for the numerical integration. 

This particular comparison showed differences of 2.2 percent of the 
temperature ratio near the corners of the passages and smaller devia- 
tions elsewhere. The temperature curves are uncertain to the same 
order of magnitude, however, because of the freedom in the extrapola- 
tion of the heat-transfer coefficient curves. The simpler relaxation 
method is therefore regarded satisfactory for the present purpose. 
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Figure 9. - Thick-walled heat exchanger with rectangular passages. 
(Shaded region is considered in calculations. ) 



Figure 10. - Thin-walled heat exchanger with rectangular 
passages. (Calculation is confined to region ABC.) 
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Figure 15. - Temperature distribution 
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